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Abstract: We investigate the infrared singularity structure of Feynman 
P> i' diagrams entering the next-to-leading-order (NLO) DGLAP kernel (non- 

singlet). We examine cancellations between diagrams for two gluon emis- 
, ■ sion contributing to NLO kernels. We observe the crucial role of color co- 

herence effects in cancellations of infra-red singularities. Numerical calcula- 
tions are explained using analytical formulas for the singular contributions. 
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Fig. 1. Two gluon Feynman diagrams contributing to non-singlet NLO kernels 



1. Introduction 

This study is part of the effort with the aim of constructing fuhy- 
exclusive (unintegrated) kernels for DGLAP [1] evolution in QCD at the 
complete NLO level. More details on this project can be found in ref. [2]. 
Let us only mention that construction of the exclusive NLO DGLAP kernels 
in ref. [2] is done following Curci-Furmanski-Petronzio (CFP) scheme [3] and 
we adopt this scheme also in our study. In short the CFP scheme uses axial 
gauge and dimensional regularization {MS) and generalizes collinear fac- 
torization developed in ref. [3]. The two particle-irreducible (2PI) evolution 
kernels Kq 

M = Co(l + Ko + Kl + ...) = CoFo (1) 

are contracted with the coefficient functions Co, which are infra-red finite. 
All infra-red collinear divergences are encapsulated in Fq, which denotes the 
sum over kernels Kq. The DGLAP NLO kernel is then extracted according 
to the scheme [3] as a single pole in Fq. More details can be found in ref. [3] 
and in ref. [2] of these proceedings. 

As it is well known individual Feynman diagrams are not gauge-invariant 
and sizeable cancellations occur among them. In particular some graphs 
contributing to NLO DGLAP kernels may contain artificial infra-red singu- 
larities, which cancel in a bigger subset of diagrams and do not appear in 
the final results. Our aim is to analyze in detail such cancellations at the 
level of the exclusive distributions, before the phase space integration, for 
the diagrams contributing to NLO DGLAP kernels in the CFP scheme. 

We shall analyze singular infra-red structure of the Feynman diagrams 
depicted in Fig. [TJ They describe emission of two gluons off a quark and 



enter into calculation of the non-singlet NLO kernel. Diagrams of Fig. 1(a) 
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1(c) are bremsstrahlung type diagrams, including interference in Fig. 1(c) 
Their C'p part is the same as in the corresponding case of QED, hence we 
shall sometimes call them "abelian" . The other diagrams of Fig. [1] include 
production of the gluon pair, see Fig. 1(f) , and its interference with the pre- 



vious bremsstrahlung diagrams. Figs. 11(d) and 1(e) Because of the presence 
of the triple-gluon vertex they are diagrams of the genuine non-abelian ori- 
gin. Moreover, since the crossed-ladder diagram of Fig 1(c) carries color 
factor equal to C|. — ^CfCa, it contributes to both "abelian" and "non- 
abelian" part of NLO kernel. 

Let us now introduce notation. The two gluon phase-space is parametrized 
using Sudakov variables: 



ki = aip + a- n + ki±, 



1,2, 



k = ki + k2, 



p — k 



(2) 



with p being the four-momentum of the incoming quark and n a light- 
cone vector. Four-vectors ki and k2 denote four-momenta of the emitted 
gluons, with their transverse parts being ki± and k2± respectively, and 
k'^ = {ki + A;2)^ being their effective mass. The sum of Qj is fixed 

ai + 02 = 1 — X. (3) 

We shall examine the distributions of two gluons in the soft limit: 




or 



ai — > 1 — X 
02 0. 



(4) 



We shall also use the "eikonal minus variables" Vi of the emitted gluons 
defined as 

The kernel is then extracted according to the scheme [3] as a single pole in 
To: 

C /a\2 f 

Fi = — y—j / d^6{l - X - ai - a2)Q{Q - max{vi,V2})p{ai,a2,vi,V2,x), 

(6) 

where C is the color factor and the function p represents contribution from 
each Feynman diagram (trace and kinematics). We shall also use the di- 
mensionless "eikonal phase-space" defined as follows: 



dai dvi , , da2 dv2 , , 
aw = d0i d(p2- 



ai f 1 



02 V2 



(7) 
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The two-gluon phase space is cut from below by means of geometrical regula- 
tor 6, namely the factors l/oj are regulated by principal value prescription: 
a" — > ^i^^2 ■ The closing of the phase space from above is ensured by the 

functior0- For the gluon pair mass we shall also use technical cut k"^ > k. In 
the numerical exercises we shall typically integrate ([6]) over the azimuthal 
angles (pi of the gluons, while concentrating on the dependence on Vi and 
aj. Also, because of the constraint in eq. ([3]), if we say that we examine the 
distribution in ai/Q2 it really means that we use ai/{l — x — ai). Also, due 
to a simple dimensional argument one of the variables Vi can be always fac- 
tored out from p function and the essential dependence of the distributions 
in variables Vi can be reduced to the dependence on the ratio y = vi/v2 
only. 



2. IR cancellations among bremsstrahlung diagrams 

In the first part we analyze "QED-like" bremsstrahlung diagrams of 



Fig. 1(a) - 1(c), We shall first show infrared cancellations among these 
diagrams in the Monte Carlo exercise and later on analyze the same cancel- 
lations analytically. 

In the following numerical exercises we keep variable x fixed and equal 
0.3. This will ensure that at least one gluon is relatively hard. 

The distributions 

f{ai,Vi)= d(l)idct>2S{l - X - ai - 02)6 {Q -ma.-K{vi,V2}) p{ai,Vi) (8) 



are plotted on Sudakov plane parametrized using variables log(ai/a2) and 
log(vi/v2), see eq. ^ for definition of Vi. 

Plots in Figs. [2] show contributions from the two ladder diagrams. They 
are obtained using Monte Carlo program FOAM [5]. As we see, their con- 
tributions appear to be strongly ordered in virtuality variables Vi of the 
emitted gluons, see for instance left part of Fig. [2j The interference diagram 
is shown in the upper right plot of the Fig. [3l It contributes in the region 
of the phase-space where both bremsstrahlung diagrams are comparable, 
which is exactly the line of equal virtualities vi = V2- The crossed-ladder 
interference diagram has a singly-logarithmic singularity along the same line 

Vl = V2- 

Fig. [3] represents both contributions: from summarized and crossed 
ladder (upper plots) together with their sum (lower plot). The singly- 
logarithmic structure, visible on the upper- left plot, disappears completely. 



^ The choice of the variable in O function closing the phase space from above can be 
different. We use max{wi,U2} or max{ai,a2} which is different from CFP choice 

= ~{p-kj-k2f. 
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Fig. 2. Two double gluon emission ladders 




Fig. 3. Summarized bremsstrahlung diagram (left), crossed-ladder diagram (right) 
and their sum (down). 



The dominant contribution is an infinite "plateau", with a long "valley" 
along the line ai — Oi2, This is not an internal infra-red singularity, however, 
but a leading-log structure. The ladder diagram is not 2PI, but quadratic 
in the leading-order kernel and therefore requires a soft counterterm to can- 
cel this doubly-logarithmic plateau, see [2]. This counterterm is necessary 
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to construct the NLO kernel and its construction is justified by the CFP 
scheme. We do not introduce it in this paper, as our aim is to present the 
universahty of gauge cancellations. Therefore we often refer to: "contribu- 
tion to the NLO kernel" rather than a kernel itself, having in mind that this 
is not a complete NLO kernel in a strict sense. Despite this dominant LO 
contribution, in the Fig. [3] one sees no structure in the soft Sudakov limit 

(SD- 

Results of the above numerical exercise can be also understood analyt- 
ically. Contributions from diagrams in Figs. 1(a) - 1(c) in the soft limit 



are proportional to simple expressions shown in Tab. [TJ Two columns in 
this table refer to two possible different Sudakov limits, with either first or 
second gluon being soft. 
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^2 (l+a;2) 2 

^F{y+xf ^ 
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p2 (1+a;') ^2..2 


(J2 (1+a;') 2 
F (y+x)'^ y 


Bx 








'^'^F(y+x)^ ^y 




SUM 


Cl (1 + x2) 


Cl (1 + X2) 



Table 1. Contributions from ladders in the Sudakov limit (up to a constant factor) 



The common denominator ^ij^^y^^i is the (rescaled) square of the virtual 

quark propagator ^ = ^(^,^^^^^^^^2 after emitting two gluons. Other factors 

come from 7-tracesE]. The last row is the sum of the three, see also the lower 
plot in Fig. [31 It is finite in the soft Sudakov limit, as the denominator 
cancels out exactly with the spinorial part of the matrix element squared. 

In the above warm-up exercise we have examined in a fine detail how 
in the NLO kernel calculations, in the axial gauge, quantum interference 
cancellations do work in practice. In fact they were exactly the same as 



^ Only "Cf-" part of the crossed ladder is taken here, hence the color coefficient is 
shown explicitly. Other unimportant factors are omitted. 
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in QED in axial gauge. Let us now turn to the genuine non-abelian gauge 
cancellations of the same kind. 



3. IR cancellations among genuine non-abelian contributions 

In this section we will discuss contributions proportional to CpC^ from 
the diagrams (c-f ) of Fig. [TJ They are of the genuine non-abelian character, 
as certified by the presence of Ca- It is interesting to see how they all 
"communicate" in the soft Sudakov limit defined below. 

We shall start with the overview of the IR cancellations in analytical 
form and next we shall illustrate them with 2-dimensional plots coming 
from Monte Carlo numerical exercises. In the following we find useful to 
use rapidity-related variables Oj in order to parametrize phase-space of two 
gluons: 

ki_L dai dai da2 da2 , 

Oj = , aW = d(pi d(p2- (9) 

Oi ai ai a2 a2 

The angular dependence enters only through the relative angle between kix_ 
and /c2± namely (/>i2, the remaining angle can be integrated out since nothing 
depends on it. 

By the soft Sudakov limit we understand that 



(10) 



while Oj is finite. The Sudakov plane in the plots will be parametrized using 
variables log(Q!i/Q2) and log(ai/a2). 

3.1. General structure of the IR non-abelian cancellations 

In Table [5] we summarize the IR cancellations of the non-abelian origin 
between the two real gluon diagrams contributions to NLO kernel in ana- 
lytical form. Formulas in Tab. [2] are the leading divergences extracted from 
two real gluon distributions in a maximally simplified form. Functions Fi 
are relatively mild and defined as follows: 

F, = ^^Fo, F2 = ^u^Fo, 
ai 02 

r— ^aia2U ^aiQ2aia2 



~ ~ 

ai(l — a2)u^ -|- 02(1 — a\) -|- 2a\a2U cos 4>\2 

where we defined u = 01/02 • Moreover we define = 1 + — 2ticos0i2, 
which is up to the 010202 factor, the effective mass of the gluon pair squared 
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Vg 
Ygl 

Yg2 
Bx 



SUM 



CfCa{1 + x^ 



CfCa{1 + x^ 



U COS 912—1^ 



2\ \ u cos (/'12 — 1 



F2 + x^Fo] 



Table 2. Cancellations of the non-abelian contributions from various diagrams due 
to gauge invariance (color coherence). Singular part for each diagram is shown in 
analytical form. Factors Fi are relatively mild, see text for their definition. 



k'^ = {ki + k2)^. Function Fq is proportional to rescaled square of the 
propagator l/g^, which is regular in the soft limit. 

The pattern of the IR cancellations in Tab. [2] is manifest. 

The most evident singularity is associated with u/a? factor, infinite when 
the effective mass of the gluon pair is zero. In logarithmic variables this u/a? 
divergence is seen as a thin infinite ridge along u = 1. It is a function of 
u and 4)12 only, strongly peaked at u = 1. In the gluonic vacuum polar- 
ization diagram this singularity is dominant. This is, however, a collinear 
singularity, not a soft one. It remains uncancelled but it is not relevant 
in our discussion of cancellation of soft singularities. We only mention it 
to understand the full singularity structure of the diagrams of interest. In 
Ygl and Yg2 diagrams the u/a? factor is present, too but the singularity is 
cancelled out by the numerators. 

The factors uFi and F2/U are functions of both u and ai/a2- They 
"soften" the sharp fall of u/o?' in the limits u — > and u ^ 00. They are 
non-zero in the soft Sudakov limit (jlOp . giving rise to a doubly- logarithmic 
IR divergenc^. These terms are present in both Vg and Yg diagrams, 
with opposite signs. This is checked and discussed below in the context of 
numerical exercises. 

In the last row of Tab. [2] the sum of all aforementioned contributions is 
presented. The terms proportional to v?F\lo? and F2I0? cancel explicitly 



The doubly-logarithmic IR divergence after phase space integration is ~ In^ j , where 
5 is cut-off variable, see Introduction. 
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Fig. 4. Triple-gluon-vertex diagrams. 

among diagram with gluonic vacuum polarization and its interference with 
bremsstrahlung. The remaining factor Fi + F2 is equal to 1 in the soft 
Sudakov limit (llOp . leaving out u/a? factor free from doubly- logarithmic 
divergences. 

Let us stress that this particular cancellation of the doubly logarithmic 
Sudakov structure of the non-abelian origin in the two-gluon distribution in 
QCD is usually referred to in the literature as the "color coherence effect", 
see for instance ref. [6]. 

The other IR divergence is caused by the presence of terms ~ Fq in 
diagrams Ygl, Yg2 and Bx, see again Table [2j Fq is nonzero along a\/a^ = 
Oii/a2-, giving rise to a single-log singularity after phase space integration. 
This has been discussed already in the case of bremsstrahlung diagrams. 
Here, however, the analytical cancellation among Ygl, Yg2 and Bx, ensured 
by ai -|- 02 = 1 — 2;, occurs in the whole phase-space, not only in the soft 
limit. 

Finally, the only singularities that remain are associated with u /o? term, 
as discussed before in this section. 

3.2. Numerical illustration of IR non-abelian cancellations 

In Fig. [Hwe show the distribution /(a^, Oi) for two gluons (averaged over 
the gluon azimuthal angles) for fixed x = 1 — ai — 02 = 0.3, from gluonic 
pair production graph (Vg) and its interference with bremsstrahlung (Yg 
= Ygl+Yg2), see also Tab.[2j Contributions from all diagrams are written 
explicitly in Tab. [21 In the plots on Figs. H]- [7] we omit their color coefficients. 
The gluonic pair production graph (Vg) has a strong peak along the line 
of equal rapiditites oi = 02 originating from u/a^ factor (off-shell gluon 
propagator). In addition, this diagram features in the plot triangular infinite 
plateau between the lines of the equal rapiditites oi = 02 and the line of 
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Fig. 5. Full cancellation of double-logs between the vacuum polarization diagram 
and bremsstrahlung-vacuum polarization interference 

equal virtualities vi = V2- It is exactly this triangular plateau which upon 
integration leads to ln^(l/(5). Note, that in Figs. S] - [7] we use variables 
log(ai/a2), contrary to the previous section, where we had log(f i/'L'2). This 
is why the line vi = V2 is now diagonal in the plots. 

Very similar, but with opposite sign, doubly-logarithmic structure is 
present in the left plot in Fig. |4] from the interference graphs Ygl-1-Yg2, 
see also Tab. [2l After adding the contributions from the above diagrams, 
see Fig. \5\ the doubly-logarithmic structure between the lines: ai = a2 and 
vi = V2 disappears. What remains, is the single-log singularity, appearing 
in a familiar shape along the diagonal line of vi = V2 (barely visible in 
the plots of Fig. H]) and collinear singularity represented here as an infinite 
ridge along the line of equal rapiditites. The latter is associated with zero 
effective mass of the gluon pair, k'^. 

The above cancellation is the well known "color coherence effect" , see for 
instance ref. |6]. It reflects the fact that the gluonic pair production graph 
(Vg) pretends that in the triangular region it is the harder gluon carrying 
octet color charge which is the "emitter", with the emission strength ~ 
CfCa- In reality, however, the emitter should be quark carrying triplet color 
charge, with the almost twice weaker emission strength ~ Cp. The role of 
the interference diagram Ygl-|-Yg2 is to correct for that and we see this to 
happen in the plot and in the formulas in Tab.O One has to remember that 
this part of the plot is already populated with the bremsstrahlung diagrams 
of the previous section proportional to ~ Cp. 

What still remains in Fig. [5] is a singly-logarithmic structure along the 
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diagonal line vi = V2- Its presence is in principle allowed in the doubly- 
logarithmic Sudakov approximation. A more subtle analysis of the soft 
limit in QCD shows that it should also vanish and this phenomenon is of- 
ten referred to as "eikonalization" , see for instance ref. [3 |8]. The job of 
bringing back the proper soft limit of the two gluon distribution and elimi- 
nating remaining single-log structure is done by the crossed bremsstrahlung 
diagram Bx (in fact its CfCa part), as it is shown in Figs. [6] and [71 see also 
Tab. El 

In right hand side of Fig. [6l the crossed-ladder diagram Bx is presented 
again. On the left hand side of Fig. [6lwe show the result of adding the triple- 
gluon vertex diagrams (Vg+Ygl+Yg2). Both plots have a characteristic 
single-log structure, seen as an infinite ridge along the diagonal line of equal 
virtualities vi = V2- However, Bx has opposite signs in its CfCa color 
coefficients. In Fig. [Tlwe see the result of adding all the above " non-abelian" 
diagrams. Bx enters with a minus sign. We see that the singly-logarithmic 
structure ~ CfCa disappears completely in the soft limit Oj 0. What 
is still present in the picture is the dominant ~ u/a^ gluon pair mass 
singularity. 




Fig. 6. The previous plot (left) and bremsstrahlung diagram (right) 



The concluding plots are shown in Fig. |S1 We presented there Feynman 
diagrams entering the kernel, both ~ Cp and ~ CfCa- In the left plot there 
are solely amplitude-squared diagrams (Brl, Br2 and Vg) and in the right 
plot all diagrams including interferences. We see explicitly the crucial role 
of "color coherence effects" (the interference diagrams) in the cancellation of 
IR singularities. In the sum of all diagrams of interest we see that remaining 
structure lies on top of the LO doubly-logarithmic plateau. The plateau does 
not enter into the NLO kernel. It is cancelled by the counterterm required 
by the kernel definition, as discussed in Section 2. 
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Fig. 7. CfCa part of the sum (Vg+Yg - Bx). 




Fig. 8. Infra-red cancellations among NLO non-singlet diagrams. Brl, Br2 and Vg 
(left), Brl, Br2, Bx, Vg, Ygl and Yg2 (right). 

4. Conclusions 

We examined the infra-red structure of the diagrams contributing to 
NLO non-singlet kernel in the unintegrated form. 

We have shown the mechanisms of gauge cancellations occurring among 
different diagrams and the importance of "color coherence effects" for this 
cancellations. 

These effects in soft Sudakov limit are examined/discussed in both an- 
alytical and numerical form. 
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